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Abstract 

The solution u{t, x, y) of the Kadomtsev-Petviashvili I (KPI) equation with 
given initial data u(0, x, y) belonging to the Schwartz space is considered. No addi- 
tional special constraints, usually considered in literature, as Jdxu{0, x,y) = are 
required to be satisfied by the initial data. The problem is completely solved in the 
framework of the spectral transform theory and it is shown that u{t,x,y) satisfies 
a special evolution version of the KPI equation and that, in general, dtu{t, x, y) has 
different left and right limits at the initial time t = 0. The conditions of the type 
J dx u{t, x,y) — 0, Jdxxuy{t,x,y) = and so on (first, second, etc. 'constraints') 
are dynamically generated by the evolution equation for t ^ 0. On the other side 
Jdxjdyu{t,x,y) with prescribed order of integrations is not necessarily equal to 
zero and gives a nontrivial integral of motion. 



1 Introduction 

We consider the Kadomtsev-Petviashvili equation ||, |^ in its version called KPI 

{ut-6uua; + uxxx)x = 3uyy, u^u{t,x,y), (1.1) 

for u{t,x,y) real. Already in 1974 ^ it has been acknowledged to be integrable since 
it can be associated to a linear spectral problem and, precisely, to the non-stationary 
Schrodinger equation 

i-idy + d^,-u{x,y))^ = 0. (1.2) 
However, to building a complete and coherent theory for the spectral transform of the 



potential u{x, y) in (1.2) that could be used to linearize the initial value problem of (1.1) 
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resulted to be unexpectedly difficult. The real breakthrough has been the discover that 
the problem is solvable via a non local Riemann-Hilbert formulation ^. 

Successively other progresses have been made. The characterization problem for the 
spectral data was solved in 0|. The extension of the spectral transform to the case of 
potential y) approaching to zero in every direction except a finite number has been 
given in ||^. The questions of the associated conditions (often called 'constraints') and 
how to choose properly in the evolution form of KPI 

ut(t, X, y) - 6u{t, X, y)ux{t, x, y) + Uxxxit, x, y) = 3d^^ Uyyit, x, y) (1.3) 

have been studied in Some additional relevant points on constraints are known, but 
not yet published |10|. 



In our opinion, the problem of the proper choice of S^T^ in (1.3) requires an additional 
investigation paying special attention to the behavior of u and spectral data at the initial 
time, say at t = 0. This study can now be done since we have at hands some additional 
theoretical tools in the theory of the spectral transform developed in |ll|, |l^ . 

In the following we will find convenient to consider the problem also in terms of the 
Fourier transform of u 

v{t,p) = v{t,p,,p2) = [ ldxdye'^^---P'y\{t,x,y), (1.4) 



(27r) 

where p denotes a 2-component vector 

P={Pl,P2}- (1-5) 



Then, the evolution equation (1.3) can be rewritten as {d q = dqidq2) 



^<*^P^ - .3^pJd\vit,q)v{t,p-q)-^Pl±^v{t,p) (1.6) 



Ot J pi 



and we are faced with the problem of properly defining the distribution Pi^ ■ The inte- 
grations here and in the following when it is not differently indicated are performed all 
along the real axis from — oo to +oo. 

According to the usual scheme of the spectral transform theory we expect that the 



spectral data evolve in time as the linear part of (1.6). This makes clear that the special 
behavior of the considered quantities at the initial time i = is just complicated by the 
nonlinearity of the evolution equation but it is, in fact, inherent to the singular character 
of its linear part. 

The final answer is that, for initial data m(0, x, y) belonging to the Schwartz space 
which are arbitrarily chosen and, therefore, not necessarily subjected to the constraint 

dxu{{),x,y) = or u(0,0,p2)=0, (1.7) 



the function u{t,x,y) reconstructed by solving the inverse spectral problem for (1.2) 
evolves in time according to the equation 

utit, X, y) - &u{t, X, y)ux{t, x, y) + Uxxx{t, x,y) = 3 dx' Uyy{t, x' , y). (1.8) 

J —too 



2 



Therefore, for smooth initial data not satisfying constraint (1.7), ut{t,x,y) has at t = 
different left and right limits and the condition 



dxu{t,x,y) = (1.9) 

is dynamically generated by the evolution equation at times t 7^ 0, and thus for these 
times we recover the result obtained in ||^. On the other side, the initial time t = Q 
requires a special investigation and the evolution form of the KPI equation we deduce in 



(1.8) is different from that proposed in ||9|. 

In the same way higher conditions ('constraints') can be obtained. But some special 
care is needed as, e.g., for the next condition we have 

dxxuy{t,x,y) = 0, i 7^ 0, (1-10) 



and the ?/-derivative cannot be extracted from the integral. In spite of all these facts 
there exists 

'dxJdyuit,x.y)=.JJdxdyuiO,x,y), (1.11) 

where the order of integrations in the l.h.s. is explicitly shown. It is clear that this 
quantity is not necessarily equal to zero and gives a nontrivial integral of motion. 
In terms of the Fourier transform v{t,p) we have 

^""^^'P^ - S^pJd\vit,q)v{t,p~q)--^Pl±^v{t^p) (1.12) 



at ^ V ^ ^ ^ ^' p^ + iot 

and 

lira v(t,pi,p2) =0, for t 7^ 0, (1.13) 
where the limit is understood in the sense of distributions in the p2 variable. 



2 The Linearized Equation 

According to the remark made in the introduction, insofar as initial value problem 

d^dtUit, X, y) + dlU{t, X, y)) = idpit, x, y), (2.1) 
U{Q,x,y)^U{x,y) (2.2) 



is concerned, we expect the linear equation and the nonlinear one (1.1) to be closely 
related. We consider here the case of initial data U[x, y) belonging to the test function 
Schwartz space S in the x and y variables. 

The study of the existence and uniqueness of the solution of the initial value problem 
defined in ( |2.lD , ( p^ ) requires to pay special attention to the behavior of the solution 
U{t, X, y) at i = 0. In fact if this function is continuously differentiable proper number 
of times in x, y and t variables and vanishing with its derivatives for (x, y) ^ 00 at any 
time, then the initial value problem has no solution if the initial data does not satisfy 
the constraint However we show in the following that if for i = (and only for 



t = 0) Ut{t,x,y) is allowed to be discontinuous and not decreasing at large {x,y) then 
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the solution [/(i, x, y) of the initial value problem exists and is uniquely determined for 
any initial data U{x, y) G S. In addition an evolution form of ( |2.l| ) satisfied by U{t, x, y) 
at all times including i = (in the sense of left/right limits) is uniquely determined. 

First of all we see that if such [/(i, a;, y) exists then, by integrating both sides of (2.1) 
for any t 7^ 0, we can rewrite (2.1) in the form 



dtU{t,x,y) = -dlU{t,x,y) + i dx" Uyyit, x" ,y), t^O 



and derive, again for any i 7^ 0, the condition 

\ dx Uyy{t, X, y) ~ 0. 



(2.3) 



(2.4) 



Let V{t,p) and V{p) denote the Fourier transforms of U{t,x,y) and U{x,y) according 
to the definition (L4). Then (2^) gives 



dt 



Pi 



and due to (2.2) 



V{Q,p)^V{p)^S. 



(2.5) 
(2.6) 



Notice that the exact meaning of p^ ^ in ( |2.5| ) is irrelevant due to condition (2.4). For 
any pi ^ the unique solution of the problem (2.5), (2.6) is given as 



V{t,p) = exp —it 



To proceed further we need the following 



Pi 



V{p). 



(2.7) 



Lemma 1 Function exp(ir/pi) defines a distribution, depending continuously on pa- 
rameter T, in the Schwartz space of the variable pi. 



i) This distribution is continuously dijferentiable in t for any t ^ : 



dr exp 







flT\ 




= — exp 






Pi 





(2.8) 



where the r.h.s. is a well defined distribution in the same space, 
a) At T = there exist right/left limits 

1 



lim — exp — 



T^±0 



Pi 



Pi 



Thus due to (2.i) we can write for arbitrary 



5^ exp — = 



Pi 



pi — iOr 



PiTiO 



IT 

exp [ — 

Pi 



(2.9) 



(2.10) 
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Proof 

i) Let us consider 

where ifipi) € S. Subtracting and adding terms we can write 

/(r) = Jdpi e"l^- [(^(pi) - I VI + V(0) / dpi e'^/^' 

Changing pi to 1 /pi in the second term we have 

/(r) = /dpi e'-/^'^ [vp(pi) - I PI |)(p(0)] + (p(O) 



(2.11) 



dpi 

|pi|>l V\ 

Now it is obvious that both terms are differentiable in r for r 7^ and 



a./(r) - z e*-/^'^ [^(pi) - t?(l- I PI 1)^(0)] + z^(0) 
J Pi 



dpi 
|pi|>iPi 



(2.12) 



(2.13) 



Substituting again pi for 1/pi in the second term we obtain that it cancels out with the 
second term in brackets. This proves i). Let us remark that ( 2.13 ) shows tha t for t 7^ 
it is not necessary to regularize the factor 1/pi in front of the exponent in (2_^). One 
can consider this factor indifferently as principal value or as (pi ± iO)~^. 
a) Again subtracting and adding terms we can write 



dpi 
Pi 



_ f dpi,r/p^ 



Pi 

sgnr (^(0) 



e'^/P^ b(pi) - i9(l- I PI 1)^(0)] + 



dpi 

\Pl\>\r\ Pi 



where now in the second term pi was substituted for rpi . To compute the limit for r — > 
notice that 

'^^^ = i^. (2.14) 



lim 



Pi 



>|r| Pi 



Thus 



lim / ^ e-/fV(pi) = / ^ [v,(pi) - ^(1- I pi 1)^(0)] ± z^(^(0) 
■^^07 Pi J Pi 



I 



— (^(pi) ±i7rv3(0), 
Pi 



where to get the second equality we chose p^^ in the integral as the principal value. Then 
( 2.10 ) follows due to (2.8) and the remark made after the proof of i). □ 

Applying this Lemma we see that V{t,p) given in (2.7), for V{p) belonging to the 
test-function space S of Schwartz, is a distribution (in the pi and p2 variables) belonging 
to the dual space S' and that V{t^p) obeys for any t the equation 



dVii.p)_ _^PA±My^t,p). 



dt 



Pi + M 



(2.15) 
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Note that the iO-term in the denominator is relevant only for t ±0 since for any t ^ 
due to rapid oscillations in (2.7) 

lim y(i,pi,p2) = (2.16) 

in the sense of distributions in p2- Thus Ut is discontinuous at i = 0. 
These properties of the distribution 

exp(^-^i^l±M^ (2.17) 

can be received as well by noting that it can be obtained as a limit in the sense of S' 
according to the following formula 



exp -it^-^—^ = lim exp -it^-^—^ . (2.18) 
V Pi J ^^+0 \ Pi+ietJ 

Then, since in the r.h.s. we have functions belonging to S for i ^ 0, we can use the 
property that the derivation with respect to i is a continuous operation in the space of 
distributions in the variables t, pi and p2 (cf. [T^ ). 

Coming back to U{t,x,y) by making the inverse Fourier transform 

U{t,x,y) = f /d25e"*«^"+^''^^-'*^'-''*^y((Z) (2.19) 



we get that it exists and solves the initial value problem (p.l|), (2.2). To study the 
properties of this solution it is convenient to extract the nonsingular part of the dispersive 
function in the exponent, i.e. the part proportional to qf. We introduce 

Uo{t,x,y)^ y"d2ge-*«i^+"J^^-*«i*F((j), (2.20) 

which solves the equation 

dtUo{t,x,y)^~dlUo{t,x,y). (2.21) 

and belongs to S in the x and y variables for any t since V{p) e S. Then from (|t|) we 
have 

Uit,x,y) ^ JJdx'dy' J{t,x - x',y- y') Uoit,x',y'), (2.22) 
where we have introduced the distribution 

J{t, X, y) = j j (fp exp ^-ixpi + iyp2 - Sit . (2.23) 

By computing the Fourier transform we get 

Jit, X, y) = ^5,(12xi - y^)-"\ (2.24) 

TT 

where the standard definition (cf. p3||) 

xl^'^=^(x)^ (2.25) 
Jx 



6 



is used. The distribution J{t, x, y) obeys the following properties 



dtJit,x,y) ^ 3dl dx' J{t,x' ,y), 



\im dt J {t,x,y) ^ ±Mi±x) dMy). 



(2.26) 
(2.27) 
(2.28) 



Now by means of ( 2.22| ) and (2.24) it is easy to check that (2. IE) gives the solution of 
(2.1) obeying all properties mentioned above. To prove that this solution is unique let 
us notice that any solution of the problem (2.5), (2.6) can differ from the one given in 
(2.7) only for a distribution localized at pi = 0, i.e. for terms proportional to S{pi) or its 
derivatives. Addition of such terms in the r.h.s. of the (2.19) violates the condition of 
decreasing of U {t, x, y) for large x. Thus in the chosen class the solution ( ^.19 ) is unique. 

Recall that we imposed to the solution to be continuously differentiable in t only for 
t 0. We have therefore to study separately the properties of ( 2.19| ) at t — 0. Due to 
(2.22), (2.21), (2.27) and (2.28) we derive that U{t,x,y) satisfies the evolution equation 



dtU{t,x,y) = ~dlU{t,x,y)+3d^y 
and that the condition 



dx'U{t,x',y), 



dxU{t,x,y) ^0, t^O, 



(2.29) 



(2.30) 



is generated dynamically. Of course (2.29) and (2.3C) are just the Fourier transforms of 
( pJ5| ) and ( ^T6|). 

In spite of (2.3C) integrating (2.22) first in y and then in x we get due to (2.23) and 



(2.24) that 



dx dyU{t,x,y) = V{0) 



dxdy U{x,y) 



(2.31) 



for any i.e. for this order of integrations the l.h.s. gives an integral of motion, which 
is not equal to zero in a generic situation. 

By means of (|^), (|^) and ( |2.24D we derive that U{t,x,y), for t {t ^ 0) and 
y fixed, decreases rapidly for x —too while for x — > too it satisfies the following 
asymptotic behavior 



U{t,x,y) 



-1 



47rv3te|a;| 
1 



dx'dy'U{x',y') 



(2.32) 



32Tr\t\V3txx'^ J J 



[ jdx'dy' [12te' + {y- y'f] U{x\ y') + o{\x\~''/^). 



This asymptotic expansion is differentiable it t and twice differentiable in y and obeys 
(2.1). Therefore we can differentiate in t the condition ( 2.30| ) and use (2.29). This 
procedure leads to the next condition 



dx X Uyy{t, x,y) = 0, t ^ 0. 



(2.33) 
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Thanks to the asymptotic behaviour of U{t,x,y) at large x derived in (2.32) one y- 
derivative can be extracted from the integral getting 



dxxUy{t,x,y) ^0, t^O, (2.34) 

but it is impossible to remove the second derivative since Jdx xU{t,x,y) is divergent. 
This procedure can be continued to get an infinite set of dynamically generated condi- 
tions. 

Equation ( p. 32 ) explains the role of constraints, i.e. conditions of the type ( 2.30| ) 



and (2.34) imposed on the initial data U{x, y), in the asymptotic behaviour of U{t, x, y) 
at large x. Precisely, for each additional constraint that is satisfied we get an additional 
x"^ factor in the decreasing law. 

Finally we notice that it is possible to weaken the conditions on the initial data. For 
example, it is enough for y(pi,p2) to be continuous at point p = Q separately in pi and 

In the nonlinear case (|l.l| ) we will consider solutions u(t, a;, y) that have the same 
mentioned smoothness properties in the x, y and t variables as U(t,x,y). The possi- 
ble discontinuity of ut{t,x,y) a,t t — can be factorized by considering the following 
representation for the Fourier transform v(t,p) of u(t,x,y) 

w(i,p) =exp (^^tt^ttM^ y(^t,p), (2.35) 

which can be considered a natural generalization of (|]^). Then v{t,p) is required to be 
continuously differentiable also a,t t — and continuous at the point p = separately in 
Pi and p2 ■ 



3 The Nonlinear Equation 

3.1 The direct problem at t = 

The direct problem has been extensively studied in ^, ^, |l^, |l^ . We report here the 
main definitions and formulae we need in the following emphasizing specific features of 
the considered problem. 

The Jost sol utio ns ^{x,y\k) are special eigenfunctions of the non-stationary Schro- 
dinger equation (1^), which are analytic in the upper and lower half plane of the complex 
spectral parameter k and whose values ^^{x,y\k) on the two sides zf^k > of the real 
k-axis 

$±(x,2/|fc) = e''''^+'''^y fi^ix,y\k) (3.1) 
are determined by the integral equation 

M±(x, y\k) = 1 + JJdx'dy' G^{x -x\y- y'\k) u{x' , y') fi^{x\ y'\k), (3.2) 

where 

1 r p-ipix+ip2y 

G^{x,y\k) = —-r d^p-— — — — (3.3) 

^ ^ (27r)2y ^ /:(p) - 2pifc T «Opi ^ ' 
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and 

C{p)=p2-pi. (3.4) 

The special role played by the spectral parameter k is stressed by separating it from 
other variables by a vertical bar. The Green function introduced in ( |3.3| ) obviously can 
be analytically continued in the upper and bottom half planes in k: 

It decreases for k ^ oo and in what follows we need the first coefficient of its l/k 
expansion. By means of (3.5) we have 

^hm kG^ix, y\k) ^ Jim ^ Jd'p _ ^ = (3.6) 

— I f ^-ipix+ip2V \ f zsz\ \ 

^ = — — +i7rsgna;%) , ±3fc > 0. 

Thus we see that the term in the expansion which is nonlocal in y explicitly depends on 
the half plane of k in which the limit is performed. Correspondingly for the first two 
terms in the l/k expansion of fj, we have 



_i_ , , , , if f f dx'dy 



- uix , y') + in dx sgn(x — x') u{x' , y) + 
J J 



Ank \ J J y — y 

°(^)' ±^^>0' (3-7) 

which also demonstrates the dependence on the half plane. 
The spectral data are defined as 

r±(a,/5) = J^jjdxdye-^^-^^-^y'^'^'--P'K{x,y)^i^{x,y\P)■ (3-8) 



Note that, due to specific dependence of the exponent in (3.8) on a and (3, spectral data 
decrease for large values of the arguments only \i a — (3 or — 0^ increase. It means 
that if we consider for some fixed parameter 7 



^ \ Iff ''■x^-iy(2f + ^ , 

P+-3^P]^Tr^ <^^dye ' V u{x,y)^i^{xMP), (3.9) 



then, since ^{x,y\(3) 1 for j3 00, 



lim ( /3 + 2 ^) = / dxdye-'^'y^u{x,y), (3.10) 



V P J (271-) 

which is not necessarily equal to zero. This remark will be crucial in the following. 

We need in the following to consider the counterpart of all these formulae in the 
Fourier transformed space. The Fourier transformed Jost solution is defined as 

„^{p\k) = jdxdye'P'^^-'P'y fi^{x,y\k), ±3fc > 0, (3.11) 
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and satisfies the integral equation 

^±(p|fc) = S^p) + !<fp'v{p~p')u^(p'\k), (3.12) 

where notation ( |3.4[ ) has been used and 5'^{p) — d{pi)6{p2)- From this equation the 
expansion in 1/fc can be easily derived. For following use we write it up to the third 
term: 

1 / £{p)v{p) 1 fj^^,v{p-p')v{p') 



(2fc)2 \ (pi ± iOp2)2 Pi ± iQp2 J p'i± iOp'2 

' ±Qk>0, k^oo. 



lim ) = w(0,27). (3.14) 



Formula (3.10) reads 

If we introduce the function 

p^{p\k)^[£{p)^2pik],.^{p\k), (3.15) 



the spectral data (3.8) can be written as 

r^{a,(3)=p^{e{a)-m\P), (3-16) 



r±(a,/3)= Jd'qv{q)i.^{e{a)~eiP)~q\p), (3.17) 
where the special two component vector 

e{k) = {k,k^) (3.18) 



has been introduced. Note that for u real (as we are considering) r'^(a,/3) = r^°'(/3, a). 

For more details the reader can refer to |jl^, where a slightly different notation has 
been used. Precisely one has to make the following identifications 

,.^{p\k) = ,.^{p,m) (3.19) 
p^{i{a) - mip) = p^Wa) - m.m)- (3.20) 

3.2 The inverse problem at t = 

The inverse problem can be solved equivalently in the (x, y) space or in the Fourier 
transformed (pi,p2) space. However, we prefer to work in the p-space because, as we will 
see in the following section, the time evolution can be more easily handled. 
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We start from the known fact that the advanced/retarded eigenfunctions i^±{p\k) of 
the non-stationary Schrodinger equation can be indifferently related to the Jost solution 
i'°'{p\k) for (T = + or for cr = — according to the formula 



iy±{p\k) = 1^" {p\k) T '2iTr J d(3 e{±a{p - k)) r-'^ (k, P) i^" [p + £{k) ~ £{(3)\(3), cr = +, 

(3.21) 

i.e. the l.h.s. is independent on the sign of cr (for a proof see g |l^, |l|, |l| ) . From (|l2l|) 
we deduce that also the sum l'-^-{p\k) + v-(j)\k) is cr independent or that 



{p\k) + ina J dp sgn{k - P) r^'^ {k, P) ly" {p + e{k) - e{P)\P) = (3.22) 
(a^-a). 

Due to this the non-local Riemann-Hilbert problem for the discontinuity of the Jost 
solution across the real k-axis (see also [|l^) is given by the equation 

iy+{p\k)~iy-{p\k) ::^-iTrY^ fdp sgn(fc - /3) r-'^(fc, /3) ly'' {p + i{k) ~ £{p)\P), (3.23) 

that has to be compared with the usual one in literature which is quadratic in the spectral 
data. 

From the analytic properties of and the Cauchy formula wc have that can 
be reconstructed from the spectral data r^{k,P) by solving the singular linear integral 
equation 

v^{p\k)=5^{p)- I I '^"^"2 r-{a,p)v''{p^l{a)- m\P) ■ 

(3.24) 

One can use the explicit fc-dependence in the r.h.s. to get a 1/A; expansion at large k 
of v^{p,k). However, before expanding the denominator in the r.h.s. it is necessary to 
perform a subtraction in the integration over a. In fact, due to ( 3.14 ) and ( 3.13 ), we 
have (for any sign of a) 

sgn(a - P) — 



lim a dp ° „ ' ' r-'^{a, P) [p + t{a) ~ e{P),P) = (3.25) 

a^oo J 2 

lim i [dp sgn.pr-'^ (a, a--] (pi + -,p2 + - ( 2a - - ) a - - 
a^ca zj \ a J \ a a \ a J a 

— sgnp2w(p), 

where to get the second line the substitution P — > a — P/a has been performed. This 
proves that a direct expan sion would furnish a divergent result. 
Therefore, we rewrite ( 3.24 ) as 

v^{p\k) = 5\p) - (3.26) 



/^^3^{E ldp'-^^^V-^).^ip^,ia)-m\P)^ 
0/ -nA sgnp2W p) \. 
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The additional term corrects the bad behavior of the integrand at large a and does not 



modify the total value of the r.h.s. of (3.24) since, once integrated over a, gives a zero 
contribution. 

Now, we can compute the first term in the expansion at large k and we get 



lim k[i^^{p\k) - S^{p)] = 

k^oo 



(3.27) 



da| Jd/3 



sgn(a - /3) — 



sgnp2w(p)}, 



±Qk > 0. 



Then, by using (3.13), we obtain 
pi ± i0p2 

- J2 sgn(a-/3) r-'^(a,/3) i^" {p + e{a) - £{P)\f3) + 



(3.28) 
sgnp2 v{p)^ 



or, taking into account that the delta distributions coming from the ±iO terms cancel 
out, we get 



Pi 



da 



{ ^ Jdp 8gn{a-(3)r-{a,(3),^^{p+e{a)-£mr3) + ^-^ sgnp2v{p)Y 

(3.29) 

where 1/pi and 1/a have the sense of principal value. 

Finally we solve the inverse problem by reconstructing the potential v{p) in terms of 
the spectral data and Jost solutions according to the following formula 



v{p) 



dapi J2 dp sgn{a- I3)r-''{a,(3)v''{p + £{a)~l{(3)\(3), (3.30) 



where we used that pi5{pi) = 0. Note that due to this fact pi cannot be extracted from 
the integral. 



3.3 The time evolution 



The time evolution in (mi) and (|3!30| ) is switched on by choosing spectral data that 
depend parametrically on the time as follows 



We need, then, to compute directly from 



(3.31) 



v^{tMk)^5^{p) 



(3.32) 



da 



a — /c =F zO 
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the evolution equation satisfied by the Jost solution {t,p\k), i.e. the second equation 
of the Lax pair. We expect it to have the form of the Fourier transform of the standard 
one 

^^.^(t,p|fc) = -2t[3{p, + k)C{p) + 2pl] v%t,p\k) - (3.33) 
2,^jj<fp'C{p')'^^v%t,p-p'\k), t^O, 
where possibly the sense of 1/p'i must be specified. Then we have to prove that 
v{t,p)^- fdapi J2 /"d/3e*'*("'-'3') sgn{a - P) r~<' {a, P) i^" {t,p + i{a) ~ l{(3)\(3) 

•' cr= + .- •' 

(3.34) 

satisfies the Fourier transform of the KPI equation. Note that for t ^ the oscillating 
term exp[4zt(a"^ — P^)] improves the behavior at large a of the integrand. Consequently, 
the potential v{t,p) satisfies the condition 

\imv{t,p)^0 for t^O (3.35) 

Pi^O 

in the sense of distribution in p2 ■ 



Taking into account the asymptotic behavior in (3.25) we see that differentiation 
of (|3|) under the sign of integral can give a divergent result at t = 0. The main 
advantage of the use of the Fourier transformed Jost solutions i/^(t,p|fc) is that they give 
an explicit possibility to remove the rapidly oscillating terms, given by the 'free' part of 
(3.33). Indeed, let us introduce the new functions 

V^{t,p\k) = e2.t[3(fc+pO£(p)+2p?] j,±(i^p|fc). (3.36) 

The kernel of the integral equation for i^^{t,p\k) does not contains a term oscillating in 
time depending on /?. In fact, by using the identity 

3{k+pi)£{p)+2pl + 2{a^~(3^) ^ 3{pi+a)C{p+£{a)-e{f3))+2{pi+a-(3f+3C{p){k-a), 

(3.37) 



we can rewrite ( p.32| ) in the form 

li^{t,p\k) = S^{p)- (3.38) 

In order to get the evolution equation for i'^{t,p\k) we have to derive this equation 
with respect to the time and exploit the exponential time dependence in the r.h.s. by 
exchanging the time derivative with the integration over a. This can be done for t ^ 
but at i = it is forbidden because of the following behavior at large a 



Sgn(Q; - P ) — 
2 

sgnp2v{0,p), 



lim a d(3 \ r-^(a, P) V{t,p + £{a) - £(/3)|/3) = (3.39) 
-5(^1) 
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that can be easily proved in analogy with (3.25) taking into account that the property 
h'^{t,p,k) S'^{p) for k ^ oo holds for any t. Note that this limit is time independent 
and determined by the initial data m(0, x, y). 
In conclusion we rewrite (3.32) as 

fda e6**-C(p)(fc-") 



(3.40) 



2{a-kTiO) 



{ J2 jd[3sgn{a~P)r-{a,(i)V^{t,p + l{a)~m\P) + 

^^(pi) sgnp2w(0,p)| 



iTT5{pi)v{Q,p) sgnp2 0{TtC{p)) 



^6itkC(p} _ 2 



k 



The last term compensates the subtracted term since 



^Q,g6it£(p)(fe-Q) 



{a =F iO) (a — fc =F iO) 
3.3.1 The first operator in the Lelx pair 



T2iTT9{TtC{p)) 



(3.41) 



The non-stationary Schrodinger equation (1.2) in the Fourier transformed space (see 
( ^ and ( |3.1lD ) reads 

[C{p)~2p,k],y{t,p\k) = [(fp'v{t,p-p')v{t,p'\k). 



(3.42) 

We want to prove, now, that the function solution of ( |3.24 ) obeys (3.42) with the 
potential v(t,p) defined in (3.34). Let us rewrite ( 3.40| ) in terms of the function 

p^{t,p\k) ^ {L{p) - 2pik] ^^{t,p\k), (3.43) 

that, due to ( 3.36 ), is related to the function p introduced in ( 3.15| ) by 

p^{t,p\k) = e^m(k+p,)cip)+2pl] p±(^t,p\k). (3.44) 

By using the identity 

C{p) - 2pik = C{p + £(a) - £(/?)) - 2(pi + a- 13)13- 2pi{k - a) (3.45) 

we get 

p^{t,p\k) = 

fda Q^itC{p){k~a) 



2(a-/fcT*0) { ^ JdPsgn{a-(3)r-{a,(3)r{t,P + e{a)-£{(3m + 

L-S{p,)\p2\v{0,p)} + 



i7r6ipi)v{0,p) \p2\e{TtC{p)) 
Pi 



„6itkC{p) _ 2 



^^g6*t£(p)(fc-a)| Idl3 sgn{a- (3)r-'^{a,(3)D"{t,p + e{a)-e{p)\P) + 

^—S{pi)sgnp2v{0,p)]. (3.46) 
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Using again ( 3.43| ) and ( |3.44 ) and the definition of tlie potential in (3.34) we have 
p^{t,p\k)=v{t,p)- 



da 



2{a-kTiO) 



Comparing this integral equation for p^(t,p\k) with the integral equation ( [3.32| ) we de- 
duce: 

p^{t,p\k) = (d?p'v(t,p~p')y^{t,p'\k) (3.47) 



that, due to (3.15), is just the Fourier transformed non-stationary Schrodinger equation 
(PI- 



3.3.2 The second operator in the Lax pair 

As we already noted, for computing the time derivative at i = of h'"{t,p\k) it is 
necessary to use the subtracted equation ( ^.4(]| ). We get different left and right limits 

^(±0,p|fc) + 



da 



dp sgn(a - /3) r-' (a,/3) ^(±0,p + l{a) - m\P) 



^ l2{a-k-iQ(j) 

C{p)[jda[ [df3sgnia-l3))7=^'{^i^^'i±0,p + £{a)~e{f3m + 



1 



a — iQ(T 



(5(pi) sgnp2w(0,p)] -2inaSipi)9{TP2<j)v{0,p)}. (3.48) 



The formula in curl brackets can be expressed in terms of v(0,p) by using ( 3.2§| ) and we 
have 



dt 

E 

-3iC{p) 



(±0,p|fc) 



da 



2{a~k- iOa) 
v{0,p) 



d/3 sgn(a - P) r-' (a, /?) ^(±0,p + ^(a) - ^(/3)|/3) 



pi ± iO 



Comparing this integral equation with (3.24) we deduce that 
d 



dt 



'(±0,p|fc) = -2i[3(pi + k)Cip) + 2pi] z/'^(±0,p|fc) - 

3i f jd^ C{p')'f^v^{±Q,p-p'\k). 
J J Pi±iO 



(3.49) 



(3.50) 



For t ^ the time derivative of v''{t,p\k) can be directly computed by using (|3.32|) . We 
get ( 3.33 ), where the distribution 1/p'i does not need to be specified as for these t the 
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potential satisfies the condition (3.35). Thus we proved that the reconstructed h''^{t,p\k) 
satisfies the evolution equation 

^yit,p\k) = -2z[3bi + k)C{p) + 2pl] iy^t,p\k) - (3.51) 



dt 



3^||d^p'C{p')^^^^'^{t,p~p'\k). 



3.3.3 The evolution version of KPI 



To derive the evolution equation for the reconstructed v{t,p) we substitute in (3.51) 
the expansion ( 3.13 ) for large k and omit terms of or der l/Zc^. Again it is necessary to 
consider separately the case t — 0. Noting in ( 3.13| ) sgn3fc = cr and multiplying the 
result by k we get 

^ ^""(±0,^)= (3.52) 



Pi + icj0p2 dt 

6ikC{p) v{0,p)\ — ^— , ^ ^ 

{pi±iv pi + im)p2 

.iL{pf +QplC(j))+^pi 



[pi + 10-0^2) 



■v{Q,p) 



3iC{p) f^,^, v{0,p-p')v{0,p') r,^, v{0,p~p') v{0,p') 



pi + ia0p2 J p'l ±iO J {p[ ± iO) {pi - p[ + ia0p2) ' 

To compensate the terms of first order in k for fc — > 00 we have to chose 

(T = ±sgnp2. (3.53) 

It means that the limits k 00 and t commute only if this condition is satisfied. 
Notice that in this case the last term in ( 3.40| ) is absent. Thus 

|(±0,p) = -zM±^,(0,p)+ (3.54) 
ot pi ± zO 



Sijd'p' [C{p)~C{p')-C{p-p 



v{Q,p-p')v{Q,p') 



p\ ± iO 
which due to the identity 

L{p) - C{p') - C{p-p') = 2p'M -Pi) (3.55) 

takes the form 

5 (±0, p) = -t^^44 ^(0, p)-6tf d^p' (pi - p[ )v{0,p-p')v{0,p'). (3.56) 



For t 7^ it is not necessary to use the subtracted form for the reconstructed v{t,p). 
Finally we deduce that v{t,p) satisfies the evolution equation 

it,P) = Ht,P) - 6« AV {Pi - p'l) v{t,p- p') v(t,p') (3.57) 



c'r ' pi+tiO 
which is the Fourier transform of the evolution equation 

ut{t, X, y) - 6u{t, X, y)ux(t, x, y) + u^xxit, x,y) = 3 j dx' Uyy(t, x' , y). (3.58) 

-too 
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